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Abstract 

We define an indefinite Wishart matrix as a matrix of the form A = W'^WT,, where E is an indefinite diagonal 
matrix and W is a matrix of independent standard normals. We focus on the case where W is L x 2 which has 
engineering applications. We obtain the distribution of the ratio of the eigenvalues of A. This distribution can 
be "folded" to give the distribution of the condition number. We calculate formulas for W real (/3 = 1), complex 
(/3 = 2), quaternionic (/3 = 4) or any ghost < /? < oo. We then corroborate our work by comparing them against 
numerical experiments. 
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Problem Statement 

Let W he an L X 2 matrix whose elements are drawn from a normal distribution. Let the two real eigenvalues of 

-*'max 



A = W WT, be denoted by Ai and A2. The condition number of A is cr 



where |Amax| = max(|Ai| , IA2I) 



and |Ajj^Jj^| = min(|Ai| , IA2I). What is the condition number distribution of A = W^WS, where S = 



Xi 



X2 



is full rank with sgn(a;2) = — sgn(a;i)? Though much is known when S is positive definite [S], to the extent of our 
knowledge, the indefinite case is rarely considered. We work out the distribution of the ratio of the eigenvalues and 
the condition number of A as it has applications in hypersensitive ground based radars [U |2] . 

A ghost /3— normal is a construction that works in many ways like a real (/3 = 1), complex (/3 = 2), and quaternionic 
{(3 = 4) standard normal [B]. In particular, its absolute value has a Xi3 distribution on [0, 00) . In general, we allow W 
to be an L X 2 matrix sampled from a /3— normal distribution but we immediately turn this into a problem involving 



real matrices namely the condition number of KSR^ , where R - 



XLI3 



X0 
X{L~1)P 



Distribution of the Ratio of the Eigenvalues 

We first write W — QR, where Q is an L x 2 orthogonal matrix and i? is an 2 x 2 upper triangular matrix |3l reduced 
QR factorization]. Note that WYW'^ = QRY.R^Q'^ and that A is similar to RY.R^ . The elements of R may be 
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chosen to be non-negative in which case it is well known that they have independent ^-distributions 



R = 



a b 






c 




X(L-l)/3 . 



and that [5] 



Xk 



x>0 



(1) 



(2) 



2fe/2-ir(fc/2)^ 

where k denotes degrees of freedom and (3 = 1 corresponds to entries being real, /3 = 2 complex, /3 = 4 quaternionic 
and general (3 for any ghost [6J. The positivity of Ri,2 merits some comment as it does not generalize well beyond two 
columns. With two columns, even for /? 7^ 1, a phase can be pulled out of the rows and columns of R and absorbed 
elsewhere without any loss of generality. 

Comment: The concept of a QR decomposition is often sensible for W^s with entries drawn from /3— normal 
distribution. As with reals, complex, and quaternions, R can be chosen to have positive diagonal elements for rank 2 
matrices, and R12 can be chosen non-negative as well, by absorbing "phases" or unit-ghosts either into Q or R'SR^ . 

This argument follows the ghost methodology described in [6 . However, mathematically the point of departure 
can be the computation of the condition number distribution for the real matrix without any mention of the 

ghosts or their interpretation. 

The joint distribution of the elements a, b,c of R with parameters L and /3 is: 



PR {a,b,c;L,l3) 



2-/3L+3„/3L-l5^-lc,3(L-l)-l 1 (^2 ^ ^2 ^ ^2)] 



(^)r(f) 



f (L-l) 



(3) 



in particular 



PR{a,b,c;L,P^l) = 
PR{a,b,c;L,l3 = 2) = 
PR{a,b,c;L,l3 = 4) = 



2a 



L-l„L-2 



exp [-i 



ttT (L-l) 

2-2L+3 (j^ _ ^) ^2L-1^^2L-3 ^_ 1 (^2 + 6^ + c^)] 

2-4(L-l) _ 1) (2i - 1) a4i-lfe3^4L^5 g^p 1 



r2 (2L) 



The first change of variables computes the matrix whose condition number we are seeking, i?Si?"^ as: R 



a^xi + b^X2 bcx2 

bcX2 



Xi 



X2 



R' = 



C^X2 





' d f 


y 





From Eq. ((T]) we see that a, 6, c are real and non-negative. The old and new 



^^7 



X2, 



f 



X2y/e/x2 



Xl 



d- 



(4) 



Note that c > implies sgn(e) = sgn(x2) and / = bcx2 sgn(/) = sgn(a;2). The Jacobian associated with this 
transformation is 



d{a.b,c) 
d{d,ej) 



4ac^xiX2 



We make the eigenvalue decomposition of the symmetric matrix 
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implying 



cos 9 — sin 9 
sin 9 cos 9 



Ai 



cos 9 sin 9 
— sin 6* cos 9 





TT' 







/ 



Ai cos^ 61 + A2 sin^ 9 
Ai sin^ 61 + A2 cos^ 9 
sm29 



The Jacobian associated with this transformation is 
ordering of the eigenvalues 



2 

d{d,ej) 



(Al - A2) 



(5) 

(6) 
(7) 

(8) 



a(Ai,A2,e) 



a;2<0<a;i => f <0 =^ Xi < < X2 
xi<0<X2 />0 ^ A2<0<Ai. 



Al — A2I. The choice of 6* £ [O, ^] nails down the 

(9) 



In summary, given 9 G [O, ^] , the constraints on d, e, f are 



a;2 < < xi =^ d > ^, e < 0, / < 0, ^ Ai < < A2 
xi <Q < X2 ^ d < ^, e > 0, / > A2 < < Al. 



Intersecting < = — ^ > with the constraint space, we obtain 



tan 



> t. 



(10) 



It is easy to see this is required by Eq. ([7]). Conversely given sgn (e) = sgn (X2) and tan^ 9 > t = —j^ we can solve for 
d and / so that d, e and / meet all of the constraints. 

The next change of variables is to write the distribution as a ratio of the eigenvalues, i.e., m = Ai, t = ^3^, and 



9. The Jacobian associated with this transformation is 

[0, 00) : X2 > 



a(Ai,A2.9) 

a{u.t,e) 



number, we integrate u over 
p{t,9;L,P) = 



[0, -00) : X2 <0 



\u\. Since we are interested in the condition 
while applying the constraint (|10l) 



2 i-xiX2)^ r (/3L) t^-^{t + if (sin 9 cos fff'^ [t cos^ 9 - sin^ 9] 



/3(L-1) 



(I) 



l(L-l) 



r(^ 



[xi (^2 cos2 6* + sin^ 9) - X2t] 



(11) 



where here and below implicitly xi and X2 are parameters as well. This in the special cases reads 

2(-xi.X2)'^/'r(L)i*-i(i + l) 



p{t,9\L,l3^ 1) 
p{t,9;L,p^2) 
pit,9;L,^ = 4) 



V^r(f)r[i(L-i)] 



t COS 



sm 



, L-l 



[xi if cos2 61 + sin^ 9) - X2t] ' 



2 al2(i-l) 



(-a;ia;2)^r(2L)t^~i (^ + 1)^ sin26l [tcos2 6'- sin 

r(i-l)r(i) [:ei (t2cos2 + sin2^?) -:E2i]''' 

(-X1X2)''' r (4L) t2^-i (t + 1)" sin^2g [t cos2 g - sin^ 9] "'"^''^ 
4r [2 (L - 1)] r (2L) [^^ (^2 eos2 + sin^ 9) - X2t\ ' 
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Lastly we integrate 6 from [0,7r/2] to obtain the distribution corresponding to the absolute value of the ratio of 
the eigenvalues p [t; L, f3). We denote the density after integrating 6 with the same symbol p without any confusion. 
In terms of the Gauss hypergeometric function [7], 2F1 (a, b; c; z) we obtain, 



p{t;L,f3) = 





)T[/3{L-1) 


+ 1] + 




l(L-l) 






\tX2 - Xif^ 



2 2 X2 — txi 



Xi — tX2 



tXl 



X2 



(12) 



In particular, 



p{t-L,P^l) 



p{t-L,p = 2) = 



■,L-1 



(-a;iX2)^/2r (i(L + 1)) r (i) {t + 1) fx,-tx 



V^r (i (L - 1)) r (i + i) \tx2-xi\' 
1 

2' 



. 2 V-^ ' 2 



tXi — X2 



2 X2 — txi 

22(L-1) _ a;2) {-XiX2 f r (L - i) (t + 1) - tX2 ' 



^ttF (L - 1) (a;i - X2) {tx2 - xi 



,2L 



tXl — X2 



-xiX2)2^r (2L + i) r (4L - 3) + 1)^ i2L-3 



0Fr [2 (i - 1)] r (4L - 1) {tx2 - xi) 

Xi — tX2 



4L 



tXi — X2 



4L 



X 2F1 4L,2;4L-1, 

\ X2- tXi 

Distribution of the Condition Number 

We can turn p{t;L,l3) to the condition number distribution by "folding" the distribution about one, i.e., take a = 
max(t, i). Therefore by folding the answer about one we can form the true condition number distribution from 
p{t;L,[3). Mathematically, 







dtp{t;L,P) = [ dtp{t;L,l3) + dtp{t;L,f3)"t:^ da (^-^p (^^; L, + p{a; L, 13)^ 
daf{a-L,(i) 

This way we arrive at the desired distribution function, f [a, L, f5), for the condition number of the indefinite 
Wishart matrices 

/ (a; L, f3) = l^p L, /s) + p (a; L, f3) , (13) 
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where p is as in Eq. (H^)- 



The condition number distribution written exphcitly is 

2/3i-i {-xiX2)^ r + i) r [/? (L - 1) + 1] {g + l)*^ a 



2 2 



/(^;i,/3) = 



where 



0iT f(L-l) r(/?L-f + l 



-a:2 

ax2 — xl 



axi — X2 



2Fi(/3L,|;/3L-| + l;r) + |^j^J^| (^)L, | ;/3L- 1 + i ) 

CT = [1, Oo) . 



In particular, 



(14) 



/(a;L,/3=l) = 



/(a;i,/3-2) = 



/(a;i:,/3 = 4) = 



2^-i(-.xi.X2)^/^r (I + ^) r (£) (g + 

0Fr(|-i)r(i + L) 

;j^|^Fi(i,L;L + i;r) + |^,;jJ^|SFi(i,L;L+i;i) 



1 (-4a;iX2)^r(L- + 
4 -X2)r(L- 1)ct4 

x2L 



o-^+^(cra;i - 2:2) cr'^+^(a;i - 0-0:2) 



{X2 - Xid)^^ {X2<7 - Xl) 



2L 



24^-1 (-a;ia;2)'^ T (2L + |) T [4 (L - 1) + 1] (a + 1)^ a 
V^r [2 {L ~ 1)] r (4i - 1) 

(^?^r=^)"2-Pl(4i^2;4L-l;r)+(,j,^)"2Fi(4L,2;4L-l;i) 



4 ^2L-3 



(15) 

(16) 
(17) 



Special Cases 

Consider the case where xi — —X2 = x which imphes r ~ —1 and the eigenvalues of A and —A are equi-distributed 
The condition number density is 



2fT^-|-ir(/3L-/3 + l)(a + l)-^(^-i) 



/3(X-l){r[| (L-l)]} 



and in particular takes the very simple form 

/(a;L = 2,/3=l) 



(18) 



7rV^(cr + 1)' 

Alternatively as /3 — cx3 the limiting condition number is non-random and equal to the condition number of the 
matrix 



y/L 1 



Xl 



X2 



Vl 1 



-1 T 



This is a simple consequence of the fact that lims-i.00 '^Xs — 1 and the condition number is independent of scaling. 
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Numerical Results 



In this section we compare the theoretical condition number distribution, Eq. (jl4p . (red curves) against Monte Carlo 
data (black dots) for the three special cases of real, complex and quaternionic matrices W. In all the plots the number 
of trials used to generate Monte Carlo data was 10^. 

Real Matrices f3 = 1 



Reals, L=2, x1=-17, x2=2 



Reals, U4, x1=7, x2=-102 




10 20 30 40 50 60 70 80 90 100 
condition number 




10 20 30 40 50 60 70 80 90 100 
condition number 



Figure 1: Real W matrices, P — 1, given by Eq. ((T 



Reals, L=6, x1=33, x2=-2 



Reals, U4, x1 =0.231, x2=-4.40 




0.03 



0.025 



0.02 



0.015 



0.01 - 



0.005 



10 20 30 40 50 60 70 80 90 100 
condition number 




10 20 30 40 50 60 70 80 90 100 
condition number 



Figure 2: Real W matrices, (3 = 1, given by Eq. (fl3|) 
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Complex Matrices /3 = 2 

Complexes, U2, x1 =-1 7, x2=2. Complexes L=4, x1 =7, x2=-^ 02 




10 20 30 40 50 60 70 80 90 100 10 20 30 40 50 60 70 80 90 100 

condition number condition number 



Figure 3: Complex W matrices, /3 = 2, given by Eq. (fT6| 



Complexes, L=6, x1 =33, x2=-2 Complexes, L=4, x1 =0.231 , x2=-4.40 
0.05, ^ ^ , , , , , 1 0-03| ' ' ' ' ' ' ' 




10 20 30 40 50 60 70 80 90 100 10 20 30 40 50 60 70 80 90 100 

condition number condition number 



Figure 4: Complex W matrices, 13 — 2, given by Eq. 
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Quaternionic Matrices /3 = 4 



Quaternions, L=2, x^=-^7, x2=2 



Quaternions, L=4, x1=7, x2=-^02 



0.05 



0.04 
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0.02 



0.01 - 




0.05 



0.04 



0.03- 



0.02 



0.01 - 



10 20 30 40 50 60 70 80 90 100 
condition number 




10 20 30 40 50 60 70 80 90 100 
condition number 



Figure 5: Quaternionic W matrices, /? = 4, given by Eq. ([T7| 



Quaternions, L=6, x1=33, x2=-2 



Quaternions, L=4, x1 =0.231, x2=-4.40 




10 20 30 40 50 60 70 80 90 100 
condition number 



10 20 30 40 50 60 70 80 90 100 
condition number 



Figure 6: Quaternionic W matrices, /3 = 4, given by Eq. ([17] 



General /3 

Here we plot Eq. for a fixed set of parameters L,xi,X2 for various ghosts. 
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L=4, x1=7, x2=-15 




Figure 7: The condition number density Eq. ((T4)) for various ghosts. As /3 cx) the distribution becomes non-random 
and can be represented by a delta function around 2.395. 



Future Work 

The natural generalization would be to extend the results to L x N indefinite matrices. We satisfied ourselves with 
N = 2 given its relevance for applications [U [5] . 
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